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ABSTRACT 

We present Very Large Array observations of ammonia (NH3) (1,1), (2,2), and CCS 
(2i - lo) emission toward the Infrared Dark Cloud (IRDC) G19.30+0.07 at ~ 22 GHz. The 
NH3 emission closely follows the 8 /nm extinction. The NH3 (1,1) and (2,2) lines provide 
diagnostics of the temperature and density structure within the IRDC, with typical rotation tem- 
peratures of ~ 10 to 20 K and NH3 column densities of ~ 10'^ cm^^. The estimated total mass 
of G19. 30+0.07 is ~1130 Mq. The cloud comprises four compact NH3 clumps of mass 30 
to 160 M0. Two coincide with 24 jjm emission, indicating heating by protostars, and show evi- 
dence of outflow in the NH3 emission. We report a water maser associated with a third clump; 
the fourth clump is apparently starless. A non-detection of 8.4 GHz emission suggests that the 
IRDC contains no bright H II regions, and places a limit on the spectral type of an embedded 
ZAMS star to early-B or later. From the NH3 emission we find G19.30+0.07 is composed of 
three distinct velocity components, or "subclouds." One velocity component contains the two 
24 /Lim sources and the starless clump, another contains the clump with the water maser, while 
the third velocity component is diffuse, with no significant high-density peaks. The spatial 
distribution of NH3 and CCS emission from G 19.30+0.07 is highly anti-correlated, with the 
NH3 predominantly in the high-density clumps, and the CCS tracing lower-density envelopes 
around those clumps. This spatial distribution is consistent with theories of evolution for chem- 
ically young low-mass cores, in which CCS has not yet been processed to other species and/or 
depleted in high-density regions. 
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1. Introduction 



Infrared dark clouds (IRDCs) are dense, cold molecular clouds seen in silhouette against 8 /urn Galactic 
plane PAH emission. Dust is relatively transparent to 8 fjm emission, so only the densest parts of molec- 
ular clouds (Ay > 50) are opaque at this wavelength dlndebetouvy et alj|2005h . Since massive stars and 
their associated clusters form in massive and dense concentrations of dust and gas, it follows that the ini- 
tial conditions of massive star formation are probably found in the highest density regions of molecular 
clouds. The GLIMPSE team imaged about two-thirds of the inner Galactic plane |Z| = 0° — 65° using the 
Spitzer Space Telescope (SST) Infra red Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8 ,um with better than 
2" res olution (IBenjamin et al.ll2003b . The GLIMPSE images contain thousands of IRDCs IPeretto & Fuller 
2009h with structure on size scales > 0.5' to as small as the IRAC resolution (see Figure [Hi. Because 
they are prominent in silhouette, IRDCs are generally expected to lie at the near kinematic distance. Pre- 
liminary st udies show that IRDCs are cold (^, 20 K) and dense (hh^ ^, 10^ cm~^), and many have masses 



> 10^ Ma, dCarev et al. 



1998 



Egan et al.l |l998] : lpillai et aPbood : bedman et allboosi : bathbome et"aDl2006 , 



20071 : IWang et al.ll2008h . iPillai et al. (2001) detected high deuteration in IRDC clumps, evidence of very low 



temperatures. However, the details of the kinematics, chemistry, and star formation within IRDCs are not 
well understood. High resolution observations of IRDCs are needed to provide information on scales small 
enough to differentiate between cloud properties in quiescent and star forming regions, a distinction crucia l 
to refining models of the earliest stages of massive star formation (|Wang et al I I2OO8I : iRathbome et al ] l2006h . 



Here we report results of a detailed study of an IRDC at radio wavelengths using the Very Large 
Array (VLA) of the National Radio Astronomy Observator>Q. IRDC G19.30-f0.07 (a(J2000)=18''25'"56\ 
5(J2000)=-12°04'30", Ylsr ~ 26 km s^^) is a filamentary IRDC with a length of ~ 2' and a width of 
~ 30". Using a flat Galactic rotation curve with Rq = 8.5 kpc and Vq = 220 km s^\ the kinematic distance 
to G19.30-i-0.07 is 2.4 kpc. We have determined the physical and chemical properties of G19. 30-1-0.07 
using ammonia (NH3) (1,1), (2,2), and dicarbon sulfide (CCS) (2i — Iq) transitions (upper level energies 
(Eu/k) of 23.8 K, 65.0 K, and 1.6 K respectively) with an angu lar resolution roughly a factor of se ven higher 
than previously reported for this IRDC, 6" compared to ~ 40" (jPillai et allboodJsakai et al.ll2008l '). To relate 
the gas properties to star formation, continuum observations at 8.4 GHz have been carried out to search for 
ionized gas associated with Hll regions in the IRDC. 

NH3 and CCS are important molecules for investigating IRDCs: both are typically considered to 
be high-density tracers (?ih2 ^ 10^^^ cm^^), and the relative stren gths of NH^ i nversi on lines can be 



used to estimate the temperature and column density of the gas (|Ho & TownesI Il983h. Further more. 
[CCSyrNH^I relative abunda nces may provide information about chemical evolution ( Suzuki et al. 19921 : 
van Dishoeck&Blakelll998k Chemical models have mainly been developed from observations of low 
mass star forming regions; this study examines whether such models might also apply in IRDCs near young 
star forming clumps. 



The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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The main goal of this paper is to examine the physical and chemical properties of the IRDC 
G 19.30+0.07 and explore the relationship between these properties and signs of star formation. The high 
spatial and spectral resolution of this study permits us to distinguish between multiple velocity components 
in the IRDC, map the temperature and density distributions in the IRDC at 6" resolution, and examine the 
relative distribution of two gas species in the IRDC. This is one of the first IRDC studies to explore the 
role of multiple cloud components, distinguished by small offsets in velocity, within an IRDC. Furthermore, 
evidence for star formation and outflows in G19. 30+0.07 enables us to examine the relationship between 
IRDC properties and star formation. 

This paper is structured as follows: in §2, we present details of our VLA observations and data reduc- 
tion. In §3, we present results of the NH3, CCS, and continuum observations. In §4 we examine the gas 
kinematics and present the derived temperature, density, and mass distributions, as well as relative abun- 
dances of NH3 and CCS in the IRDC. In §5 we discuss the relationship between physical and chemical 
properties, kinematics, and star formation within the IRDC. Conclusions are given in §6. 

2. Observations and Calibration 

Table [T] summarizes the parameters of the VLA observations. NH3 (1,1) and (2,2) and CCS (2i — Iq) 
emission at ~ 22 GHz were observed with the VLA in D configuration. The 3. 125 MHz bandwidth contained 
the main and innermost satellite hyperfine lines of the NH3 (1,1) and (2,2) transitions. Two pointings were 
used to fully cover the highest opacity regions, as determined from the 8 /urn images. All data were calibrated 
and imaged using AIPS. The CCS observations were conducted during the EVLA upgrade, and eight of the 
twenty-seven antennas had been converted to EVLA antennas. Data were obtained using both the VLA and 
EVLA antennas. Doppler tracking could not be used during the upgrade period, so the AIPS task CVEL 
was applied during calibration to correct for motion of the Earth relative to the local standard of rest during 
the observations. In the 22 GHz datasets, the uncertainty of the absolute flux densities is estimated to be 
10%. We estimate the absolute positional accuracy to be better than 1". 

Images were made of the two pointings at each frequency using natural weighting of the u-v data and 
CLEAN deconvolution. The AIPS task FLATN was used to mosaic the two pointings and apply a primary 
beam correction. The NH3 (1,1) and (2,2) mosaics were cut off at the 50% power level of the primary beam. 
The synthesized beam in the full-resolution NH3 images is ~ 4"; however, the NH3 data were later smoothed 
to 6" resolution to increase the signal-to-noise ratio prior to fitting the emission lines with Gaussian profiles 
as a function of velocity (see §4). The multiplication factors needed to convert from Jy beam^' to beam- 
averaged brightness temperature (K) for the various images are listed in Table [T] 

Because the CCS emission is weak and extended relative to that of NH3, a 30 kX Gaussian u-v taper was 
applied while imaging the data, resulting in a resolution of 7"7 x 5"7. The u-v taper was chosen to maintain 
the highest resolution possible while still providing a significant (> 3a) CCS detection. Additionally, a 
more conservative primary beam cutoff at the 65% power point was used while mosaicing the CCS so that 
noise peaks at the edges of the primary beam were not amplified into apparently significant detections at the 
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low contour levels. The 65% power point cutoff was determined by reducing the primary beam cutoff until 
spurious 3a detections no longer appeared at the image edges. 

Continuum emission at 8.4 GHz was observed with the VLA in A configuration. These data were also 
calibrated and imaged with the AIPS reduction package, using the model for the flux calibrator supplied 
with the package. The uncertainty of the absolute flux density scale is estimated to be 10%. Robust u-v 
weighting and CLEAN deconvolution were used during imaging. The resolution of the resulting image is 
~ 0"3. 



_We have also re-reduced and analyzed the G19.30+0.07 VLA H2O maser data published by' Wang et al. 



(|2006i) (project AW659; the full details of these observations are described by Wang et al.). The data for the 



two pointings toward G19. 30+0.07 were reduced in the usual manner using AIPS. After self-calibration, 
the RMS noise in a single 0.66 km s^' wide channel is 20 mJy beam^^ toward the NE pointing and 
15 mJy beam^' toward the SW pointing These noise level s are about five times better than those reported 



for the ensemble of sources observed by IWang et al.l (12006b . The synthesized beam for these maser data is 
approximately 1"5 x 0"9. 



3. Results 



3.1. NH3 



The NH3 (1,1) main Une emission, integrated over 24. 1 to 30.3 km s ^ is shown overlaid on an infrared 
3-color Spitzer image of G19.30h-0.07 in Fig. [I] The three c olors in Fig. [Dare GLIMPSE 4.5 ^im (blue) and 
8 jum (green) emission and MIPSGAL dCarey et all l2009h 24 /urn emission (red). Fig. [U shows that the 
NH3 (1,1) emission is spatially well correlated with the 8 fum extinction. NH3 (2,2) emission was also 
detected throughout the IRDC, although due to the fainter line emission the (2,2) emission was not as 
well defined near the IRDC boundaries. The RMS noise in the NH3 (1,1) integrated intensity image is 
5 mJy beam^^ km s^^ 



In our discussions of IRDC structure, we use the guidelines described by iBergin & Tafallal (|2007h : 
clouds have masses of ~ 10^ to 10^ Mq and radii of ~ 1 to 10 pc, clumps have masses of ~ 10^ to 10^ M© 
and radii of 0.1 to 1 pc, and cores have masses of ~ 1 to 10^ M0 and radii of ~ 0.01 to 0.1 pc. Clumps 
produce star clusters, and cores produce individual stars. A 20" FWHM Gaussian unsharp mask was used 
to identify the compact structure in G 19. 30-1-0.07. The selection criteria for NH3 clumps were peaks having 
structure < 20" with integrated intensities > 25 Jy beam^^ km s^' in the masked image. G19.30h-0.07 
contains four apparent NH3 clumps, labeled in Fig. [T] as CI, C2, C3, and C4. CI and C2 coincide with 
24jm emission and have previously been studied at mm wavelengths with ~ 1 1" resolution (|Rathbome et al 
20061) . C3 is not associated with any 24 ^m emission; however, a H2O maser is located nearby (see §3.4). 
C4 is an NH3 emission peak not associated with any IR emission or masers. The 24 ^m source S I (Fig. [B is 
not co incident with a peak in the NH3 emission. There is a H2O maser with velocity 19 km s~^ (|Wang et al. 
20060 spatially coincident with SI. This maser's velocity is slightly offset from the velocity of G19.30+0.07, 
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~ 26 km s^^ While it is possible that the maser is associated with an unrelated object along the line of sight, 
the probability of two dense cores capable of forming a H2O maser along the line of sight to SI is small. It 
is more probable that the velocity offse t between the ma ser and the IRDC is caused by motion within the 



IRDC related to star formation activity dWang et al.ll2006h . 



Multiple velocity components and broad line-wing structures are evident in the NH3 emission, as shown 
in Figure |2] Double-line profiles indicate locations in the cloud where there are two overlapping kinematic 
components. The double-line profile features cannot be attributed to self-absorption: as demonstrated in 
Fig.|2j), the double-line profile structure is seen in the more optically thin (2,2) emission as well as the (1,1) 
emission. Furthermore, the double-line profiles are not correlated with NH3 high optical depth (see §4.2). 
Finally, the individual components of the double-line profiles can be traced continuously from regions with 
single-line profiles at the same velocity; it is in overlap regions that they appear double, as illustrated by 
the profiles shown in Fig. |2] There are three distinct velocity components over the entire G 19. 30-1-0.07 dark 
cloud, which we will henceforth refer to as "subclouds." The spatial extents of these subclouds are derived 
from fits to the NH3 line profiles, as described in §4.1. The subcloud positions relative to the 8 /um emission 
are shown in Fig. |2] 



3.2. CCS 

Where both CCS and NH3 emission is detected their velocities match very well. However, in many 
locations there is emission from one molecule and not the other. Figure |3]shows the distribution of the CCS 
and the NH3 (1,1) main line emission along with representative CCS profiles; peaks in the CCS emission 
have been labeled as P1-P5. The CCS in Fig.[3]is integrated over four channels (1.2 km s^^) centered on the 
velocities of the subclouds as determined from the NH3 emission. The NH3 emission in Figure |3] has been 
smoothed to match the resolution of the CCS emission {I'.'l x S'.'l) and integrated over the same velocity 
ranges. The RMS noise in the CCS integrated intensity images is 1.3 mJy beam^^ km s^^ 



3.3. 8.4 GHz Continuum 



No 8.4 GHz continuum emission is detected from G19. 30-1-0.07. The RMS noise in the field is 
0.025 mJy beam^^ A 5a (0.125 mJy) upper limit on any free-free continuurn emission imp lies an up- 
per limit on the UV photon flux of A/^c ^ 8 x 10"^^ s~^ ( Mezger & Henderson 1967 ; Condon|[l992r) . However, 
it is unlikely that all of the emitted UV photons go into ionization of H and He. The UV flux could be 
reduced by absorption by an accretion disk or dust surrounding the protostar. Thus, the UV flux limit is 
likely an underestimate of Nc- Furthermore, early-stage H II regions can be small, dense, and optically thick, 
and consequently have low radio continuum flux densities. Further discussion of the protostellar mass limit 
implied by the UV flux Umit is given in §5.3. 
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3.4. H2O Masers 



From our re-reduction and analysis of the Kvang et al.l (120061) H2O maser data we have detected two 
regions of maser emission. A previously unreported maser with a peak flux density of 1.58 Jy was detected at 
a(J2000)=18''25'"51.95^ 5(J2000)=-12°05'14"0 at an LSR velocity of 28.8 km s^i. This maser is located 
within the C3 ammonia clump (Fig. [T]), and is the f irst reported evide nce of star formation associated with 
this clump. The second H2O maser, first reported in IWang et al.l (|2006|) . has a peak flux density and velocity 
of 0.51 Jy and 19.0km s-\ respectively, at aposition of a(J2000)=18''25'"58.64\ 6(J2000)=-12°04'2i;'5. 
This maser is located close to the 24 /urn source SI. The reported peak flux densities have been corrected for 
primary beam attenuation. 



4. Physical Properties of the Dense Gas 
4.1. Gas Kinematics 

The line profiles from the spatially smoothed 6" x 6" resolution NH3 (1,1) and (2,2) images have been 
fitted using a non-linear least squares Gaussian fitting routine, providing central velocities, linewidths, and 
amplitudes. Fitting errors for the NH3 central velocity and linewidths are typically < 0.1 km s^^ and 
< 10% for the amplitude, with the uncertainty reduced by approximately a factor of two in areas with a high 
signal-to-noise ratio. 

The (1,1) line fitting was done in two steps. First, the main (1,1) Une was fitted. Second, the main 
line and two satellite hyperfine lines were fit simultaneously with a triple-gaussian profile, using the results 
from the first (main-line only) fit to set initial guesses for, and constraints on, the fits to the satellite lines. 
The initial guesses and constraints were set as follows: (1) the initial guesses for the central velocities of 
the satellite lines were set at the velocity of the main line ±7.7 km s^\ constrained to ±0.6 km/s of this 
initial guess; (2) the initial guess for the width of the satellite lines was set to the width of the main line, 
constrained to ±0.6 km s^^ The fit to the (2,2) line was not constrained by the results of the (1,1) fit. 

In areas with double-line profiles (see Fig.|2l), the fitting included additional steps. One of the (1,1) main 
lines and its associated satellite lines were fitted as described above. Then the second (1,1) main and satellite 
lines were fitted as described above. Finally, the results of these fits were used to set the initial guesses 
and constraints for a simultaneous (six-gaussian) fit that separated the various emission contributions. The 
double-profile fits of the satellite lines used the same main-line based constraints as described for single- 
profile fits. The fits to (2,2) double-line profiles were done simultaneously using a two-gaussian profile; as 
in the single-profile case, the fits were not constrained by the (1,1) line fits. CI and C2 exhibit an emission 
component with broad line- wings, possibly indicating outflows; in these areas the line- wing component was 
modeled as an additional Gaussian during simultaneous line fitting. 

Three subclouds were identified based on fits of the NH3 (1,1) profiles. Their median LSR velocities 
are centered at 25.7 km s^\ 26.7 km s^\ and 28.4 km s^K The spatial distributions of the velocities and 
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linewidths for each of the subclouds are shown in Figure H] Internal motion is indicated within each subcloud 
by central velocity changes of < 1 km s^^ shown in Figs. |4^, |4j), and|4]:. Typical FWHM linewidths are 
~ 1.2 km s"^ but the linewidths are ~ 0.6 km s^' at the subcloud boundaries and increase to > 2 km s^^ 
near C 1 and C4, shown in Figs. |4}i, |4^, and|4f . Since C4 contains no evidence of star formation, the increased 
velocity dispersion at C4 may be due to a nearby outflow originating from C2. 

As discussed in §3.2, CCS emission is found in all of the subclouds. The CCS emission lines could 
not be fitted reliably because of their low signal-to-noise ratio, but based on the profiles shown in Fig. |3]we 
estimate the CCS FWHM to be ~ 0.6 km s^^ and therefore assume the CCS emission from each subcloud is 
contained within 1.2 km s^' centered at that subcloud's velocity. Where both CCS and NH3 are detected, 
the hnewidths and emission line structure from both molecules are in good agreement, even in areas with 
multiple velocity components. Thus, CCS velocity components were separated by averaging over 1.2 km s^^ 
centered at the subcloud velocities determined from NH3. 

NH3 (2,2) line profiles were used to examine the kinematics of the line-wing features, as the greater 
separation between the main line and satellite hyperfine components for the (2,2) transition simplified si- 
multaneous fitting of narrow and broad components. The line- wings toward CI and C2 had typical FWHM 
linewidths of ~ 12 km s^^ A residual image cube was made of the NH3 emission with the fitted narrow 
component subtracted, so that only the broad component remained. The integrated blue- and red-shifted 
emission was obtained from the residual data by integrating over the velocity ranges 14.6 to 26.6 km s^^ 
and 26.6 to 38.7 km s^^ Representative profiles and the red- and blue-shifted integrated NH3 (2,2) emission 
are shown in Figure |5] 

The red and blue line- wing emission from CI is spatially coincident with, and centered on, a 24 /urn 
source, and can be interpreted as an outflow oriented close to the line of sight. The infall velocity derived 
from the mass and radius of CI (see §4.4) is ~ 2 km s^^ insufficient to account for the broadening observed 
in the line-wing emission. The age of Cl's outflow cannot be determined from the flow velocities because 
of the line-of-sight orientation. The line-wing emission at C2 can be interpreted as an outflow originating 
at the 24 /urn source, with the red shifted outflow extending ~ 0.4 pc NE and the blue shifted outflow 
extending ~ 0.3 pc SW of the source. The outflow from C2 appears well collimated, with a collimation 
factor (length/width) > 5. The outflow's inclination angle is unknown, so we assume an angle of 45° in 
estimating the measured length of the outflow. The central velocity derived from a Gaussian fit to the broad 
component at each pixel within the C2 outflow gives an average outflow velocity of ~ 5 km s^^ corrected 
for an assumed inclination angle of 45°. The outflow's kinematic age was determined by dividing the extent 
of the outflow by the average outflow velocity, resulting in an age of ~ 10^ yr. 



4.2. Optical Depth, Temperature, and Density Structure 



NH3 inversion transitions have hyperfine splitting due to magnetic dipole i nteractions and the in terac- 
tion of the nitrogen electric quadrupole moment and the electron electric field (IHo & Towneslll983h . The 
optical depth of the NH3 (1,1) main line, X(i 1^), can be derived from the ratio of the main and satellite 
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hyperfine line strengths. Assuming the excitation temperature is the same for the main and satelhte lines, 
the ratio of the brightness temperatures j,„,) and can be expressed as 

Th(l.l.m) _ 1 -exp[-T(i.i.„)] 



where 'm' and '5' indicate the main and satellite hyperfine components, Tj, is the measured brightness tem- 
perature, X(i 1 is the optical depth of the main quadrupole hyperfine component of the (1,1) line, and a 
is the theoretical intensity ratio of the satellite line to the main line if the hyperfine levels ar e populated ac 
cordi ng to their statistical weights. For the innermost satellite hyperfine transition, a = 0.28 (IRydbeck et al. 

mi. 



The total optical depth of the (1,1) line, X(i 1), can be determined from '^(i.i,m) ^rid the relative intensities 
of the hyperfine components, such that X(i i) = ^X(i j ,„), where ^ is the ratio of the total (1,1) line optical 
depth to the optical depth of the main hyper fine component. If th e NH3 (1,1) hyperfine levels are populated 
according to their statistical weights, C, = 2 (IRydbeck et al.ll 19771 ). 



The optical depth X(i 1) ranges from ~ 1 to 4 at the boundaries of the subclouds and increases towards 
the cloud centers. Near CI, C2, and C3 X(i j) is ~ 5 to 10. The uncertainty in X(i 1) is < 30% where the 
NH3 signal to noise is highest. At the subcloud boundaries, the uncertainty in X(i jj increases, and is on 
the order of a factor of two. At positions where X(i j) > 20, 7^(1, i^v) approaches r^(i j,„,) too closely to yield 
accurate results; at these positions a X(i i) value of 20 was assumed. The optical depth distribution is shown 
in Figure |6] It is notable that areas where NH3 double-line profiles are detected, shown as subcloud overlap 
in Fig. |2l are not coincident with the highest optical depth gas, ruling out self-absorption as a cause of the 
double peaked line profile. 

The ratio of the brightness temperatures of the NH3 (l,l,m) and (2,2,m) lines can be used to determine 
the rotation temperature Tk(\,\ : 2,2) that characterizes the population distribution among the (1,1) and 
(2,2) states (Iho & Towneslll983l') : 



7)^(1,1:2,2) 





0.284 , 




In 




'^(l.l.m) 



^b(2.2.,m) 
Tb(\.\.m) 



(1 -exp[-T(i^l_„,)]) 



(2) 



The distribution of Tr for locations where the statistical uncertainty in Tr is 10% or less is shown in Figure |7] 
The gradient is steep near the NH3 clumps, and consequently Tr, in particular, suffers from the effects 
of beam dilution caused by smoothing the data. Thus, we created temperature maps at the original ~ 4" 
resolution of CI and C2, where the NH3 emission is strong enough to obtain good hne fits. The higher 
resolution images are shown outlined in black in Fig. Ha). 

Values of Tr range from 10 K up to ~ 30 K, above which r^(2,2,m) approaches 7),(i 1 „j) too closely 
to yield reliable temperature estimates. A value of 35 K was assumed at points where r/,(2,2,m) ~ ^fe(i,i,m) 
to indicate points with uncertain temperatures. The higher resolution images of the clumps show that Tr 
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reaches this maximum at CI and SI, and 7); ~ 20 K in C2, indicating heat sources at these locations. 
Heating is consistent with the presence of 24 fxm emission in CI, C2, and SI. An area with > 35 K 
approximately 10" east of CI has no corresponding peak in NH3, and the nature of this peak in Tr remains 
unknown. 



Using earlier results from Danby et all (11988|) . 
tween IRDC gas kinetic temperature, Tk, and Tr, 



Tafalla et al. 



(|2004f) empirically find a relationship be- 



Tk 



Tr 



l-§ln 



l-lexp(Ti^) 



(3) 



This relationship is good to better than 5% when Tr is 5 to 20 K. Above 20 K, the relationship is not as 
accurate. Since the majority of the gas in G19. 30+0.07 has Tr < 20 K, we apply this relationship and find a 
mean Tk in the IRDC of ~ 17 K compared to a mean T/; of ~ 15 K. 

The column density of the NH3 (1,1) transition was calculated from the total optical depth and bright- 
ness temperature of the (1,1) transition, X(i 1) and 7;,(i.i), which include all hyperfine components. In the 
optically thin limit, j.m)' where ^ = 2 and is defined as above. The column density of the NH3 

(1,1) transition was determined using 



N, 



AV 



(1,1) 



(1,1) 



y8E(2) (l-exp[-T, 



(1,1) 



(4) 



where AV is the FWHM linewidth (in velocity units), v is the frequency of the (1,1) transition, and A(i 
is the Eins tein coefficient for the sum of all hyperfine transitions. Equation 4 contains an optical depth 
correction (IGoldsmith & Langerlll999h that corrects for the opacity of the (1,1) emission. The total column 
density is obtained using the NH3 partition function, Q[Tr]: 



N(li)Q[Tr] 
Ntotal = exp 

^1,1 



1,1 



klR 



(5) 



Equation 5 assum es the same Tr among all NH3 metastable levels, with negligible contribution from non- 
metastable levels feachiller et al. 1987 ). £'(1,1) is the energy above the ground state and g(i.i) is the total 
statistical weight of the (1,1) inversion state. Values for th e partition function were derived from JPL line 
catalog values calculated at 9.375 K, 18.75 K, and 37.50 K jpickett et al.ll 19981) using a linear interpolation. 
Equations 4 and 5 assume that the emission fills the synthesized beam, and neglect the microwave back- 
ground. In areas where the excitation temperatures approach the cosmic microwave b ackground temperatur e 
(2.7 K), Equation 5 may underestimate the column density by roughly a factor of two (iBachiller et al. Ill987b . 
Corrections for the opacity of the (1,1) emission were included in the calculation of the NH3 column densi- 
ties; however, the column densities will be lower limits in areas of very high optical depth (x(l, 1) > 20), 
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since the optical depth cannot be measured accurately in these areas. Fortunately, such high-opacity regions 
cover a very small area of the IRDC, and are not associated with the clumps. 



The NH3 column density distribution is shown in Figure [8] The column densities lie in the range 10^^ 
to 1.6 X 10'^ cm^^. In all three subclouds, there appears to be a lower density envelope (~ 10^^ cm^^) 
around higher density clumps (-^ 4 x lO'^ to 1.6 x lO'^ cm^^). The uncertainty in column density is < 30% 
within the clumps, but at the cloud boundaries the values are only good to approximately a factor of two. 
CI, C2, and C3 are coincident with column density peaks. 

Pillai et all \200d ) found the [NH3]/[H2] abundance ratio in IRDCs ranged fr om 0.7 x 10"*^ 



10 X 10"^. We assume a [NHi l/rH?] abundance r atio of 3 x 10'^ (IHariu et al.l ll993b. The H2 column 



to 



densities imply peak Ay > 100 (|Bohhn et al.lll978l : ISchultz & Wiemedll975h . in good agreement with the 
observed 8 /urn absorption. To convert from column density to volume density, the line of sight depth must 
be determined. We assume a cylindrical geometry for the filamentary IRDC, with the semi-major axis 
("length") directed from the northeast comer to the southwest comer of the IRDC. We approximate the 
semi-minor axis ("width") of the IRDC, aligned perpendicular to the semi-major axis, to be ~ 0.35 pc. The 
assumption of cylindrical geometry was then used to calculate the depth at each point in the IRDC. The 
calculated H2 number densities then rang e from 3 x 10^ to 5 x 10^ .9"^"^; values consistent with previ 
ous n umber densities reported for IRDCs ( Carey et al. 1998 : Egan et aDll998 : Pillai et al. 2006; Wang et al. 
2008k 



4.3. CCS Column Density 

The CCS (2i — lo) column density, N\, is determined using 



8%^\^kgi 1 AV 
Ni= , , --^Th- 



hc^ g2A2i '^81n(2)' 



(6) 



where the symbols are defined as in Equations 4 and 5. Equation 6 assumes the CCS (2i — Iq) emission 
is optically thin and fills the synthesized beam. We assume ThAV to be the integrated emission of the CCS 
over 1.2 km s^^ Using the partition function Q[Tex\, where T^x is the temperature describing the relative 
level populations, the total column density is obtained: 



total 



- exp 



(7) 



where E is the energy of the Iq level above the ground state. Val ues for the partition function were interpo- 
lated from JPL hne catalog values in the same manner as for NH3 ( Pickett et"aDll998 ). We assume a CCS T^^ 
of 10 K based on typical values of NH3 Tr in areas where CCS was detected. Typical CCS column densities 
ranged from lO''' to 3 x lO'^ cm^^, roughly two orders of magnitude lower than the NH3 column densities. 
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Unfortunately, the unknown CCS Tg^ introduces significant uncertainty in the column density measurements; 
using an excitation temperature of 5 K reduces the column densities to 5 x 10^^ to 10^^ cm~^, while 15 K 
produces column densities ranging from 2 X lO^^ to 6 x 10" cm . Due to the large uncertainties in the 
CCS excitation temperature and linewidth, these CCS column densities are likely accurate to only an order 
of magnitude. 



The CCS T„ can be constrained using CCS (43-32) observations of iSakai et all (I2008h (HPBW ~ 37"). 
Sakai et al.l (120081) used a non-detection of CCS (43-32) emission to report an upper Umit on the total CCS 
column density of < 3.7 x 10^^ cm^^ at CI and < 3.5 x 10'^ cm^^ at C2. We used these results to derive 
an uppe r limit for the CCS (43-32) column density, at CI and C2. To compare our observations with 
those of ISakai et al.l (120081) we smoothed our CCS A^i column density image to m atch the Nobeyama Radio 
Observatory resolution, FWHM 37", and measured N\ at the ISakai et al.l (|2008h pointings. By comparing 
our measured N\ values and the A'^3 upper limits and assuming local thermodynamic equilibrium, we derived 
a CCS Tf.^ upper limit of < 5 K. 



4.4. Mass 



IRDC and clump masses were derived from NH3 column densities. The mass of the IRDC was deter- 
mined by integrating the NH3 column density (shown in Fig. [D over the entire source area (clumps as well 
as diffuse gas), giving the total number of NH3 molecules in the IRDC. The number of molecules was con- 
verted to total ma ss by assuming the mass is dominated by H2 and using the assumed [NH3]/[H2] abundance 
ratio of 3 X 10"^ dHarju et al.lll993h an d an H? mass of 3. 3 4 x 10'^"^ g. The IRDC mass is ~ 1 130 Mq, con- 
sistent with the ~ 900 Mq reported by IPillai et all (|2006l) . bathbome et al.l (|2006 l) determined a total IRDC 
mass of ~ 400 M© using 1 .2 mm dust emission. The apparent discrepancy between the masses determined 
from NH3 observations and dust observations may be attributed to the large systematic uncertainties in the 
dust emissivity, dust to gas ratio, and [NH3]/[H2] ratio. 

The individual masses of CI, C2, C3, and C4 were calculated by integrating NH3 column density over 
only the clumps' spatial extents, and assuming an [NH3]/[H2] abundance ratio of 3 x 10^^. The clump 
radii were determined by fitting the NH3 emission distribution in CI, C2, C3 and C4 with a 2-D Gaussian 
using AIPS task JMFIT and using the FWHMs as the radii of the clumps. Clumps CI, C2, C3, and C4 
have masses 140, 160, 30, and 30 M©, respectively. The mass uncertainties given in Table |2]are based on 
uncertainty in the column densities and clump radii; it should be noted that uncertainties in the distance to 
the IRDC and the [NH3]/[H2] a bundance ratio introduc e additional systematic uncertainty. While the mass 
uncertainty may be quite large. iRathbome et al.l (l2()06h de termined masses of 113 for CI and 114 M0 
for C2 using 1 .2 mm dust emission. 



Rathbome et al 



(120061) assumed a dust temperature of 15 K, comparable 
to the tempera tures we have measured in this IRDC (mean Tk ~ 17 K, mean Tr ~ 15 K). The clump radii 
determined by IRathbome et al.l (|2006|) were 0. 15 pc for CI and 0.21 p c for C2, comparable to the ~ 0.1 pc 
determined for C 1 and C2 in this study. The IRathbome et al.l (|200a) masses are i n good agreement with 
those found in this study, despite the fact that the overall IRDC mass determined by IRathbome et all (|2006r) 
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is a factor of three lower than that determined in this study. The fact that the masses agree on the scale of 
the clumps but not on the scale of the cloud may be due to sensitivity and/or resolution differences between 
the two studies. 



5. Discussion 
5.1. Clump stability 

We examined the stability of C 1 , C2, C3, and C4 by comparing their gas mass to their virial masses. The 
virial mass (Mv,>) is the mass at which thermal plus turbulent gas pressure balances gravitational forces in the 
clumps, i.e., the gravitational and kinetic energy are in equipartition. The ratio of virial mass and clump mass 
{Myir/Mciump) is a basic indicator of whether the clumps are gravitationally bound and likely to collapse: if 
Myir/M clump < 1. the clump contains enough mass to collapse under gravity. This balance analysis assumes 
negligible contributions to clump stability from additional forces such as magnetic fields and clump rotation, 
neither o f which are known for th e clumps in G19.30+0.07. Additionally, clump confinement by external 
pressure dBertoldi & McKeelll992h has not been included. We calculated M,,„- from the clump radius {Rdump, 



determined in §4.4) and the NH3 FWHM hnewidth (AV) using 



The NH3 (1,1) FWHM linewidth was measured at the density peak of the clumps. Values of AV, M^ir, and 
Mvir/M clump are listed in Table |2] 

The clumps all have Myir/Mdump > 1> inconsistent with indications of star formation that imply the 
clumps are undergoing gravitational collapse. However, Equation 8 assumes virial equilibrium. If the 
clumps have already collapsed and harbor outflows the conditions of virial equilibrium are violated, and 
the linewidths will not accurately predict the thermal and turbulent motion supporting the clumps. Indeed, 
clumps CI, C2, and C3 show evidence of outflow through either broad NH3 linewidths or H2O masers. 
Clump C4 has Myir/Mdump ^ 8; this high ratio may be a result of increased linewidths toward C4 caused by 
the outflow originating from C2 (see §4.1). It is also worth noting that where the cores are optically thick the 
AV may be overestimated, which would lead to a greater overestimate of virial mass (which is proportional 
to AV^) than of clump mass (which is proportional to AV). 



The virial masse s of CI and C2 rep orted by IPillai et all (l2006h are 893 and 823 Mq, significantly larger 



than our M^r values. 



Pillai et al 



(120061) used larger clump radii and higher linewidths in their virial mass 
calculations because of the lower resolution of their study (40"). Our higher resolution shows most of the 
mass is contained within a smaller radius than previously thought. Furthermore, we w ere able to r esolve 
velocity gradients in the emission near the clumps (Fig. The larger beam used by IPillai et all (l2006h 
would have confused this emission, contributing to their larger linewidth measurements. 
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5.2. The Relationship Between NH3 and CCS 



Chemical evolution models comparing CCS and NH3 distributions have mainly been developed 
from observations of low mass star forming regions and the environs of H II regio ns; we now examin e 
whether such models rn ight a lso apply in IRDCs. Chemica l models such as those of ISuzuki et alj (119921) . 
Bergin & Langeri (Il997h . and Ivan Dishoeck & Blakel (Il998h predict that CCS traces chemically young gas 
in lower density, quiescent envelopes of molecular clouds, while NH3 traces dense, evolved clumps. Time- 
dependent gas-phase modeling suggests that for low-mass star-forming clumps the [CCS]/[NH3] abundance 
ratio is a measure of the age of clumps, i.e., a "chemical clock": CCS abundance peaks in chemically 
young gas and decr eases as the CCS compo nents are reprocessed into molecules with stronger binding en- 
ergies, such as CO. Millar & Herbsti (| 199(1) and iNejad & WagenblastI (| 19990 propose CCS is reprocessed 
on timescales of ~ 10^ years, while de Gregorio-Monsalvo et al. (2005, 2006) concluded that CCS has a 
lifetime of ~ 10^ years after the onset of star formation. 

Freeze-out onto dust grains likely plays a significant role in chemical changes during clump evolu- 
tion. Dust grain chemistry is particularly important for explaining the high abund ances of NH3 observed 
in dense clumps (^ Bergin & Langer 1997 : Bergin & Tafalla boOT : Sakai et al. 2008 ). When CO freezes out 
onto dust grains the abii ndance of NH3 is enhanced, because gas-phase CO destroys NH3 parent molecules 
(|Bergin & Langeiill997h . Furthermore, the abundance of hydrogenated m olecules, including NH3, may be 
enhanced as ices evaporate off dust grains dvan Dishoeck & Blake 1998 ) near hot clumps in star-forming 
regions. Freeze-out of CCS onto dust grains ma y also be responsible for low CCS abundances in gas clumps 
where ?ih2 > 10^ cm~^ (IBergin & Lange Dust grain chemistry also explains how the chemical clock 

may be reset. iDickens et al. (2001) propose collisions in gas with properties similar- to IRDC subclouds can 
trigger MHD waves and/or grain impacts that provide sufficient energy for thermal desorption of molecules 
from grain surfaces, leading to the presence of early-time molecules like CCS. 

The critical densities of the NH3 (1, 1) and CCS (2i — In) emission lines differ by almost two orders 



of ma gnitude, ~ 2 x 10^ cm"^ for NH3 (IDanbv et al.lll988r) and ~ 10^ cm"^ for CCS (IWoUcovitch et al. 



19971) . which can complicate the interpretation of molecular distribution. However, despite its higher critical 
density, we detect CCS emission in the lower density envelopes around the high-density NH3 clumps rather 
than toward the center of the dense clumps, suggesting their spatial distributions are the result of chemical 
evolution in the IRDC rather than a density effect. 

Previous studies of the che mistry in G 19. 30-1-0.07 have been limited by poor resolution. Prior to this 
investigation, IPillai et al.l (120061) observed G19. 30-1-0.07 in NH3 emission using the Effelsberg 100 m tele 
scope with 40" resolution, but lacked the resolution needed to comment on small-scale structure. ISakai et al. 
J2008h observed N2H+(l-0), HC3N(5^), CCS(43-32), CH30H(7-6), and NH3 (1,1), (2,2), and (3,3) emis- 
sion towa rds the clumps in G 19.30-i-0.07 with the NRO 45-m telescope (HPBW ~ 37" for CCS and ~ 73" 
for NH3). ISakai et al.l (120081) report abundances suggesting the clumps in G19. 30-1-0.07 are more chemically 
evolved than low-mass starless cores. The high spatial resolution data we present here improve upon these 
findings by showing how morphology, temperature, and density structure in the IRDC impact the distribu- 
tions of the NH3 and CCS. 
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Fig. [3]shows the integrated CCS emission in each of the subclouds and compares it to the NH3 inte- 
grated emission. Although CCS emission is associated with the subclouds at 25.7 km s^^, 26.7 km s^^ and 
28.4 km s^^ the spatial distribution of the CCS and NH3 emission are not coincident. Enhancements in 
the CCS with clump-like morphology are labeled as PI to P5 in Fig.O no evidence of NH3 clumps exist at 
these locations. 

\i3 



Typical CCS column densities ranged from 10 cm to 3 x 10 cm 



>13 tn 1 10^^ f-m--i \Ye cxamlncd the 

[CCS]/[NH3] column density ratio in areas with detectable emission from at least one of the molecules. 
In CCS peaks lacking detectable NH3 emission we assign an NH3 column density upper limit of 10^^ cm^^. 
Where NH3 but not CCS was detected we assign a CCS column density upper limit of 10^^ cm^^. At clumps 
CI, C2, C3 , and C 4 the [CCS]/[NH3] ratio is < 10"^ slightly higher than the upper limit de t ermined by 



Sakai et al.l (120081), and consistent with active, star forming clumps (IBergin & Tafallal 12007: ISakai et al 



20081 



h. At CCS emission peaks, [CCS]/[NH3] ratios are > 0.02, consistent with results from Isuzuki et al. 



(I1992I ). who found [CCS]/[NH3] > 0.03 in low-mass starless clumps. A map of the column density ratios is 
shown in Figure |9] 

Comparing Fig. |9]with Fig. |3] shows that PI is associated with the [CCS]/[NH3] peak and located 
in the subcloud located at 25.7 km s^\ the only subcloud lacking evidence of star formation. The CCS 
morphology and [CCS]/[NH3] ratio in this subclou d are consistent with a chernically young clump in a 
stage too early to present evidence of star formation (|de Gregorio-Monsalvo et al.ll2006r) . 



The CCS peaks P2, P4, and P5 are located at subcloud interfaces, and P3 is spatially coincident with the 
outflow from C2. The spatial coincidence of these CCS peaks with subcloud-subcloud and outflow-subcloud 
interfaces indicates collisions may be resetting the chemical clock of the gas at these locations, consistent 
with the iDickens et all (1200 ih model. 



One major caveat must be considered in the analysis of CCS emission. CCS (2i — Iq) is a low-lying 
energy state transition (upper energy level only 1.6 K above ground), so it is difficult to kn ow whether 



enhanced CCS emission is due to excitation effects or increased abundance. Using the results of ISakai et al. 



(120081 ) we were able to roughly constrain the CCS excitation temperature to < 5 K. However, the CCS 
(43 — 32) observations were low resolution and only provided a column density upper limit. Additional 
transitions of CCS must be measured at comparable resolution to this study to more accurately determine 
the CCS excitation temperature. Sensitive interferometric observations at millimeter wavelengths could 
explore whether the abundance measurements in the ground state are consistent with abundances measured 
from higher energy transitions. 

We draw the following conclusions from the CCS detections: (1) CCS and NH3 are present in each 
of the three subclouds comprising G19. 30-1-0.07, but each molecule traces different parts of the subclouds; 
(2) CCS appears to trace the boundaries of the IRDC while NH3 peaks in the inner parts of the IRDC, 
consistent with CCS enhancement in low-density envelopes and NH3 enhancement in denser clumps, in 
spite of the fact that the CCS transition has the higher critical density. This emphasizes that NH3 does 
not give the full picture of the dense gas in star-forming regions; (3) in some cases, CCS appears to have 
clump-like morphology. The [CCS]/[NH3] peak at PI is located in the only subcloud lacking evidence 
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of star formation. This supports the hypothesis that CCS cores may represent the very earliest stages of 
core formation; (4) CCS peaks at subcloud-subcloud or outflow-subcloud interfaces may indicate collisions 
are resetting the chemical clock in some parts of the IRDC; and (5) further investigations at higher-energy 
transitions are needed to distinguish excitation and abundance effects in CCS emission. 



5.3. Star Formation in G19.30+0.07 



Although early studies of G 19.30+0.07 did not find any evidence of star formation (iCarey et al.lll998h . 
recent investigations, including this one, indicate substantial active star formation in this IRDC. The ob- 
served evidence for star formation is: (1) archival MIPSGAL 24 ^m observations showing areas of hot dust 
associated with CI, C2, and SI; (2) structure in the NH3 spectral line observations consistent with outflows 
originating from CI and C2; (3) H2O masers coincident with C3 and SI (Fig. [Hi. 

No 8.4 GHz continuum emission was detected in the IRDC. A 5a detection limit of 0.125 mJy beam 



sets the upper limit o f the UV photon flux at < 8 x 10^^ s ^ , implying a stellar type of B 3 or later dSmith et al. 



2OO2I ; |Panagialll973h . Howev er, if the protostar is accreting material, much of the stellar UV radiation will 
be absorbed close to the star (IChurchweIllll997h . Also, in an IRDC it is likely that some UV photons are 
absorbed by dust. Thus, stars with types hotter than B3 could be present in G19.30-f -0.07 but be undetecte d 
by our measurements. In fact, Orion Sourc e I, an early-B star w ith mass > 8 M0 (IMatthews et al.ll201( 
has a n 8.4 GHz flux density of ~ 1.1 mJy (jMenten & Reidlll995r ) at a distance of ~ 400 pc (|Menten et al. 



20071). If Orion Source I were located at the distance of G19.30-i-0.07, 2.4 pc, its flux density would be only 
~ 0.03 mJy, below the detection limit of this study. 

Clump masses of > 10^ M0 an d the presenc e of outflows suggest that CI and C2 may be protoclusters. 
Using the initial mass function from lScald (119861) and assuming a 50% star forming efficiency, these clumps 
would form stars up to masses of '-^^ 3 Mq. The smaller NH3 clump and H2O maser without accompanying 
24 ^m emission at C3 suggest C3 likely contains protostars that are either lower-mass or younger than those 
in CI and C2. The 24 fjm emission toward SI suggests the presence of at least one protostar, consistent 
with the H2O maser coincident with SI. It may be SI is more evolved than CI and C2, and has destroyed 
or dispersed any dense gas associated with its formation. However, no shell of dispersed NH3 is detected 
around S 1 . 

Assuming the linewing structure discussed in §4.2 indicates outflows, the outflow masses provide in- 
sight into the masses of the stars contained in C 1 and C2. The column densities in the outflows were deter- 
mined as in Equations 6 and 7, using the integrated emission from the NH3 (2,2) line-wings and assuming a 
Tex of 10 K for NH3. The column densities were summed over the outflow area to determine the masses of 
the outflows. The red and blue outflow components of CI each contain ~ 2 M0 and the components in C2 
each contain '-^ 1.6 Mq, suggesting intermediate or high mass protostar progenitors. However, we note that 
NH3 only traces high density gas, so observations of the low-density gas are needed to better constrain the 
total outflow masses. 



Various models attempt to explain the formation of massive stars, but recent work suggests accretion 
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in turbulent clumps is the likely creation mechanism (IBeuther et aLll2007h . iBeuther et all (120071) present an 



evolutionary sequence of massive stars starting from high mass starless clumps: low- and intermediate-mass 
protostars form within the massive starless clumps and accrete matter until becoming a massive protostellar 
object, which develops into a massive star and clears the natal gas from its environment. C 1 and C2 may be 
examples of clumps containing the embedded accreting protostars proposed in this model. 

Based on the results from this single IRDC, the role that subclouds play in massive protostar and 
protocluster development is still unclear. They may indicate an early stage of fragmentation, cloud forma- 
tion/collision, or a superposition along the line of sight. The subclouds appear distinct, but if less dense gas 
below our detection limits connects the subclouds, they may all be part of a larger structure. Further work is 
needed to understand the role of subclouds in IRDCs. 



5.4. Is G19.30+0.07 a "Typical" IRDC? 



Understanding the role of IRDCs in star formation requires a survey of their properties. We have 
examined one IRDC in great detail, but is G19.30-i-0.07 a typical IRDC, and how general are our findings? 
Our results are consis tent with previous studies that indicate G19.30-h0.07 has temperatures and de nsities 
typical for an IRDC dcarev et al.lll998l : Ipillai et al.ll2006l : Iwang et allboosl : iRathbome et"aDl2006h . The 
high resolution of this study compared to previous work has revealed new IRDC features. Notably, the 
kinematics have been analyzed in detail. We find that the IRDC is composed of three subclouds, each of 
which has a slightly different velocity. The universality of multiple subclouds in IRDCs and the role they 
play in protocluster evolution requires further investigation. 



6. Conclusions 

Massive star formation continues to be poorly understood, largely because early stages are difficult to 
observe due to the short timescales involved and the embedded nature of massive protostars. Because of their 
large masses, low temperatures, and high densities, IRDCs are currently the most likely candidates to host 
the first stages of massive star formation. Understanding their properties may shed light on the processes 
that lead to the creation of massive protostars. Although IRDCs are identified using IR observations, their 
high opacity makes radio wavelengths an important probe of their structure. We have examined IRDC 
G19.30-i-0.07 using VLA spectral hne observations of the NH3 (1,1), (2,2), and CCS (2i - Iq) transitions. 
From this work, we draw the following conclusions: 

1. G19.30-i-0.07 has structure down to the resolution limit of our observations (4" to 6"). At a distance 
of 2.4 kpc, the structure is < 1.4 x lO'^ AU in size. This level of detail is seen in the IR images as well 
as the NH3 and CCS observations. 

2. G19. 30-1-0.07 is composed of three subclouds at different velocities. Star formation and [CCS]/[NH3] 
levels may indicate that one of the subclouds is at an earlier evolutionary stage than the other two. 
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Further investigation is needed to determine if IRDCs are typically composed of multiple subclouds, 
and the role subclouds play in IRDC evolution. 

3. Typical values of Tr in G19.30+0.07 range from 10 to 30 K with a mean value of 15 ± 3 K. Increases 
in at CI, C2, and SI indicate heat sources at these locations. 

4. NH3 column densities in G19. 30+0.07 are ~ 10^^ to 1.6 x 10^^ cm^^, with lower density envelopes 
around higher density clumps. 

5. The NHa-derived masses of CI, C2, C3, and C4 are 140 ±20, 160 ±40, 30 ±20, and 30 ±20 Mq. 
The clumps have masses significantly lower than their virial masses; however, signs of star formation 
indicate CI, C2, and C3 are gravitationally bound and collapsing. In clumps that are collapsing or 
contain outflows, M^triai /Mcore is not a reliable metric for testing gravitational boundedness. 

6. CCS and NH3 are present in each of the three subclouds comprising G19. 30+0.07, and appear to 
trace different parts of the subclouds. The NH3 emission closely traces the IR extinction seen in the 
GLIMPSE 8 /7m images, and peaks where clumps are forming and evolving in the IRDC. Conversely, 
the CCS tends to trace the boundaries of the IRDC, where NH3 c olumn densities indicate the IRDC 



is more diffuse. Given the high critical density of CCS (~ 10^) (IWoUcovitch et al.lll997h . the CCS 
detections at the IRDC boundaries demonstrate that NH3 does not give the full picture of the dense 
gas in star-forming regions. 

7. In some places, the CCS appears to have clump-like morphology. The CCS peaks may show the very 
earliest stages of clump formation and/or result from colhsions that reset the chemical clock of the 
gas. 

8. Stars are forming in G19.30+0.07, most obviously in clumps CI and C2. Younger or lower mass 
clumps may look like C3, where a H2O maser but no IR emission is present, or C4, where no signs 
of star formation have yet been detected. SI, where a H2O maser coincides with a 24 /an source but 
no NH3 clump is detected, may be a more evolved protostar, or associated with NH3 at a velocity 
outside our observed bandwidth. We did not detect 8.4 GHz emission from the IRDC, although the 
8.4 GHz emission from the early-B star Orion Source I would be below our detection limits if it were 
at a distance of 2.4 kpc. 
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Table 1 . VLA Observation Details 



Parameter NH3 (1,1) NH3 (2,2) CCS (2i - Iq) 8.4 GHz 



Observation Date 


2005 Dec 10, 11 


2005 Dec 10, 11 


2007 Feb 19 


2006 Jan 30 


NRAO Proposal ID 


ADS 16 


AD516 


AD556 


AD524 


Rest Frequency (GHz) 


23.6944955 


23.7226333 


22.3440330 


8.4351 


VLA Configuration 


D 


D 


D 


A 


Total Bandwidth (MHz) 


3.125 


3.125 


3.125 


50 


Number of channels 


128 


128 


128 


1 


Spectral Resolution (kHz/km s^^) 


24.4/0.3 


24.4/0.3 


24.4/0.3 




Primary Beam FWHM 


1.9' 


1.9' 


1.9' 


5.3' 


Synthesized Beam 


4:'5x3"4 


4:'6x3:'4 


7"7x5"7 


0"31xO"19 


RMS Noise (mJy beam^' channel^ ^) 


4.5 


4.5 


3 


0.003 


Beam Position Angle 


4.8° 


-0.6° 


6.7° 


20.1° 


Conversion Factor (K/Jy) 


144 


138 


55.7 


2.9 xlO^ 


Smoothed Synthesized Beam* 


6"0x6;'0 


6"0x6f'0 






Smoothed Beam Position Angle* 


0° 


0° 






RMS Noise (Smoothed Data, mJy beam^^ channel^') 


6 


6 






Smoothed Conversion Factor (K/Jy) 


60.5 


60.3 






Largest Detectable Angular Scale 


94" 


94" 


100" 


16" 


Phase and Amplitude Calibrator 


J1832-105 


J1832-105 


J1832-105 


J1832-105 


Flux Density Calibrator 


3C48 


3C48 


3C286 


3C48 


Bandpass Calibrator 


3C273 


3C273 


3C345 





*A11 NH3 data used in line-fitting and calculations have been spatially smoothed to 6" resolution to increase the signal to noise ratio. The 
only un-smoothed NH3 data are presented in Fig. [T]and in the high-resolution portions of Fig. |7^. It should be further noted that in Fig. |3l the 
NH3 data have been smoothed to match the CCS beam, 7"7x5"7. 
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Table 2. NH3 Clump Properties 





CI 


C2 


C3 


C4 


Mclump (H2* ) (Mo) 


140 ±20 


160 ±40 


30 ±20 


30 ±20 


Mclump (dust**) (Me) 


113 


114 








0.11 ±0.002 


0.12 ±0.003 


0.07 ±0.006 


0.09 ±0.005 


AV (km s-i) 


2.8±0.1 


1.8±0.1 


1.6±0.1 


2.1 ±0.2 


Mvirial (M©) 


540 ±50 


260 ±40 


110±20 


240 ±40 


^^virial/^clump 


4±1 


2±1 


4±2 


8±5 



*Derived from NH3 as discussed in §4.4. 



*Dust results are quoted from Rathbome et all (120061) . 
~ 0.012 pc assuming a distance of 2.4 kpc. 



-22- 



04:00.0 



o 
o 
o 

CM 



o 



u 



30.0 - 



Q-12:05:00.0 



30.0 - 




18:26:00.0 58.0 56.0 54.0 

Right ascension (J2000) 



52.0 



25:50.0 



Fig. 1. — 4.5 /rni (blue), 8.0 ^um (green), and 24 /jm (red) 3-color image of IRDC G19. 30+0.07. White contours represent the 
integrated intensity from the main hyperfine component of the NH3 (1,1) line. The NH3 (1,1) resolution is 4"5 x 3"4 and RMS 
noise at the center of the primary beam is '--^ 5 mjy beam^' km s^' . Contours mark 4o, 80, 12<J, 16CJ, and 20CJ and then increase 
in 8CT intervals, where C is defined at the center of the primary beam. Blue lines mark the 50% power point of the primary beam; 
noise level s at the edge of th e field are twice those at the center of the primary beam. The light blue asterisks indicate H2O maser 
detections 1 Wang et^boOfii . this paper). The four NH3 clumps identified by this work are labeled CI, C2, C3, and C4. 
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Fig. 2. — Three representative line profiles of NH3 (1,1) and (2,2) emission from three different positions in G19. 30+0.07. The 
solid and dot-dashed lines show the (1,1) and (2,2) data, respectively. Profile (a) shows a single velocity line profile, observed 
at a(J2000)=18''25"'52.2', 5(J2000)=-12°05'20!'7. Profile (b) shows emission from multiple velocity components, separated by 
~ 1.5 km s-1, observed at a(J2000)=18''25'"52.P, 5(J2000)=-12°04'59!'9. Profile (c) is from a point in the northern ammonia 
clump, a(J2000)=18''25'"58.8', 5(J2000)=- 12°04'02"9, showing a broad line-wing component. The locations of the representative 
profiles are indicated on the GLIMPSE 8 /jm image. The line profiles were used to identify three subclouds, separated by velocity 
offsets ~ 1 km s^'. The subcloud locations are shown overlaid on the GLIMPSE 8 /jm image; the subcloud velocities are centered 
at 28.4 km s"' (red), 26.7 km s"' (yellow), and 25.7 km s"' (blue). 
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Fig. 3. — CCS (2i - lo) distribution in IRDC G19.30+0.07 and representative profiles of CCS emission. The greyscale is the 
8.0 /jm image from GLIMPSE. Red contours represent CCS (2i — Iq) emission integrated over different velocity components: (a) 
25.2 to 26.5 km s"', (b) 26.2 to 27.4 km s"', and (c) 27.7 to 28.9 km s"'. The angular resolution of the CCS image is 7"7 x 5;'7. 
The CCS RMS noise is 1.3 mjy beam^' km s^' . CCS increments start at 3a and increase in 2(T increments, where a is defined at 
the center of the primary beam. The white contours show where the detection power of the primary beam in the CCS observations 
drops to 65%. Blue contours represent NH3 (1,1) emission smoothed to match the CCS resolution {I'.'l xS'.'l) and averaged over 
the same velocity ranges. The NH3 RMS noise is 1.8 mJy beam^^ km s^^, and the contours indicate emission levels of lOo, 30CJ, 
60cj, and lOOo. Arrows indicate CCS peaks, labeled PI, P2, P3, P4, and P5. Representative profiles toward PI, P2, and P4 are 
shown in (d), (e), and (f), respectively. 
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Fig. 4. — Velocities and FWHM linewidths within tlie three subclouds that make up G19. 30+0.07. Pixel spacing is 3", resolution 
is 6". Black outlines indicate the boundaries of the different subclouds identified in Fig.|2] The LSR velocities are shown in (a), (b) 
and (c), and the NH3 (l,l,m) FWHM linewidths are shown in (d), (e), and (f). Velocity components are centered at 25.7 km s^' in 
(a) and (d), 26.7 km s"' in (b) and (e), and 28.4 km s"' in (c) and (f). 
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Fig. 5. — Broad line-wings in addition to narrow NH3 emission are found near CI and C2. Three representative profiles of the 
NH3 (2,2) emission are shown from different locations in G19. 30+0.07. In (a), (b), and (c), the total (2,2) line emission is shown 
as a dashed line and the residual emission (line-subtracted, as described in the text) is shown as a solid line. In (a) and (b), broad 
line-wing components are present; for comparison, (c) lacks line-wing emission. The locations of the representative profiles are 
indicated on the GLIMPSE 8 /jm image (d). Contours in (d) show integrated intensity from the NH3 (2,2) residual images. The 
blue and red-shifted line-wing emission was integrated over 14.6 to 26.6 km s^' and 26.6 to 38.7 km s^^; these velocity ranges are 
indicated by red and blue arrows in (a), (b), and (c). Blue contours represent the integrated intensity blue-ward of the NH3 (2,2) line 
center, and red contours represent the integrated intensity red- ward of the NH3 (2,2) line center. The la uncertainty in the emission 
integrated over 6 km s^' is ~ 0.025 Jy beam^' km s^' . Contours show 2o, 6ct, and then increase in 3CT intervals. 
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Fig. 6. — The optical depth of NH3 (1,1) in the three subclouds is shown. The components are centered at 28.4 km s 
26.7 km s^' (yellow), and 25.7 km s^' (blue). Contours start at t(1,1)=3 and increment by 3. 



(red). 
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Fig. 7. — The NH3 rotation temperature distribution in the three subclouds is shown. Only points with uncertainty < 10% have 
been mapped. Grey outlines indicate the boundaries of the different subclouds. The component centered at 25.7 km s^' is shown 
in (a), 26.7 km s^' in (b), and 28.4 km s^' in (c). The resolution in the images is 6" except in the areas outlined in black in (b), 
which have a higher resolution of 4'.'5 x 3'.'4 and pixel spacing O'.'S. Arrows in (b) indicate the locations of CI, C2, and Si. 
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Fig. 8. — The NH3 column density integrated over all velocities is shown in (a). The column density contribution of the separate 
velocity components is shown in (b). The separated velocity components are centered at 25.7 km s^' (blue), 26.7 km s^' (yellow), 
and 28.4 km s^' (red). Contours indicate lO'^ cm^^, 2 x lO'^ cm^^, and then increase in steps of 2 x lO'^ cm^^. 
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Fig. 9. — Column density ratio map showing the relative distribution of CCS and NH3 ([CCS]/[NH3]). Only areas with detectable 
emission from at least one of the molecules are mapped. In CCS peaks lacking detectable NH3 emission we assign an NH3 column 
density upper limit of 10^^ cm^^. Where NH3 but no CCS was detected we assign a CCS column density upper limit of lO'-' cm^^. 
Contours scale logarithmically, with log([CCS]/[NH3]) = —2.9 (violet), —2.7 (cyan), —2.5 (green), —2.3 (yellow), —1.9 (red), and 
-1.7 (white). 



